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Abstract 

Using molecular dynamics computer simulations, we investigate the dynamics 
of the rotational degrees of freedom in a supercooled system composed of rigid, 
diatomic molecules. The interaction between the molecules is given by the 
sum of interaction-site potentials of the Lennard-Jones type. In agreement 
with mode-coupling theory (MCT), we find that the relaxation times of the 
orientational time correlation functions C^i\t), C^\t) and Ci(t) show at low 
temperatures a power-law with the same critical temperature Tc, and which 
is also identical to the critical temperature for the translational degrees of 
freedom. In contrast to MCT we find, however, that for these correlators the 
time-temperature superposition principle does not hold well and that also the 
critical exponent 7 depends on the correlator. We also study the temperature 
dependence of the rotational diffusion constant and demonstrate that at 
high temperatures is proportional to the translational diffusion constant 
D and that when the system starts to become supercooled the former shows 
an Arrhenius behavior whereas the latter exhibits a power-law dependence. 
We discuss the origin for the difference in the temperature dependence of D 
(or the relaxation times of C^^^) and D^- Finally we present results which 
show that at low temperatures 180° flips of the molecule are an important 
component of the relaxation dynamics for the orientational degrees of freedom. 

PACS numbers: 61.43.Fs, 61.20. Ja, 02.70.Ns, 64.70.Pf 
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I. INTRODUCTION 



In the last ten years mode coupling theory (MCT) [|T]-|^ has led to a strong interest 
in the phenomenon of the glass transition, since this theory makes detailed predictions on 
the dynamics of glass formers in the vicinity of the glass transition and thus challenges 
the experimentalists to test these predictions. Starting from the microscopic equations of 
motion, MCT derives, by using certain approximations which are believed to be quite reliable 
for simple liquids, an equation of motion for the density correlator. One of the main results of 
MCT is the existence of a dynamical transition at a temperature T^, at which, in the so-called 
idealized version of the theory, the dynamics of the system undergoes a transition from an 
ergodic (T > Tc) to a nonergodic behavior (T < Tc), and which can be interpreted as a glass 
transition. Due to the presence of ergodicity-restoring processes, commonly called "hopping 
processes", most glass formers show to a certain extent deviations from the predictions of 
the idealized theory, since this version of MCT does not take into account these sort of 
processes. In the extended version of the theory hopping processes are taken care of and 
one finds that in the vicinity of a crossover in the behavior of the dynamics can still be 
observed Many of the predictions of MCT were confirmed by experiments and in 

numerical simulations on various glass forming systems. For reviews the reader is referred 
to Refs. [0^. 

In the real world most of the good glass formers are molecular systems. One of the 
important differences between molecular systems and simple liquids is that the former have 
orientational degrees of freedom (ODOF). The dynamics of these ODOF can be measured, 
e.g., by dielectric spectroscopy, light scattering, and NMR. For the molecular glass formers 
salol, e.g., it has been shown that MCT gives a satisfactory description of the relaxation be- 
havior of these systems 0. This conclusion seems to contradict the outcome of experiments 
in which the dielectric response of this system was probed [0 because of the following fact: 
MCT predicts that for all observables that couple to the density fluctuations the imaginary 
part of the corresponding susceptibility exhibits a minimum at the same frequency Wmin, if 
one is very close to Tc. Light scattering experiments have shown that in salol this minimum 
occurs at a frequency of about 2GHz if T = 253 K ||^. However, Dixon et al. showed that 
no such minimum is observed in the dielectric function e"{ijj) for frequencies up to 18GHz 
for T = 255 K [0 thus seemingly contradicting the conclusions of Ref. that MCT gives 
a correct description of the dynamics of this system. Similar conclusions were drawn for 
glycerol. However, in recent extensive dielectric measurements on glycerol this minimum 
was found and it was shown that the data can be described by the /5-correlator of MCT 
for simple liquids although in this experiment the location of the minimum is indeed 
different from the one in the light scattering experiment. This could be due to at least 
two reasons: First, the MCT-equations for molecular liquids might be different from those 
for simple liquids, hence leading to different predictions. Second, the predictions of MCT 
for simple liquids can be used, but corrections to the asymptotic laws of MCT have to be 
taken into account. These corrections make the reduction theorem invalid and therefore the 
location of the minimum depends on the observable |T0[. We also note that very recently 
evidence has been given that a minimum exists in e"{uj) for salol also [0. However, its 
location remains still undetermined 

The above discussion shows that the role of the ODOF for the glassy dynamics in the 
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supercooled regime is certainly not settled and thus remains an interesting field of research. 
One possibility to gain some insight into the dynamics of supercooled liquids is to perform 
molecular dynamics computer simulations of such systems |]12[ . Simulations of simple liquids 
have shown that this method can indeed be very useful to understand the dynamics down 
to temperatures close to Tc [|T^-[T5|. For supercooled molecular liquids, however, only few 
numerical simulation have been done so far. Signorini et al. investigated the dynamics of 
the ionic glass former [Ca{NO'i)2\oA[KNOz\f).Q (CKN), where the NOz unit was treated as 
a rigid molecule and Sindzingre and Klein studied methanol ||T^ . OTP was investigated 
by Lewis and Wahnstrom |T^,T^, who modelled the molecule with a rigid isoscale triangle of 
Lennard- Jones particles, and by Kudchadkar and Wiest |2y] who used a 18-site, three-ring 
model. Very recently Sciortino et al. presented their results on supercooled H2O (also taken 
as a rigid molecule) [pT] , p2| . In these papers a two-step relaxation process, as predicted 
by MCT, was observed for both, the translational and rotational degrees of freedom. But 
a more quantitative analysis in the framework of MCT was essentially restricted to the 
translational degrees of freedom (TDOF). In particular it does not become clear from these 
studies whether the ODOF freeze at the same temperature as the TDOF do since no detailed 
analysis with respect to this point was made (although very recently evidence was given, that 
for water the freezing temperatures of the ODOF and the TDOF are very close together [0). 
The investigation of this point is one of the major objectives of the present work. For this 
we will focus on the time scale of the a-relaxation and discuss our findings regarding the 
/^-relaxation elsewhere 
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We also mention that in a very recent paper Ma and Lai 
investigated the dynamics of the translational degrees of freedom for a dumb-bell shaped 
molecule and argued that the molecular character leads to a decrease of the MCT parameter 
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The molecular system we consider consists of molecules in which two different atoms 
are connected rigidly. Apart from diatomic molecules with head-tail symmetry, this is the 
simplest molecular system one can choose. Although the choice of a linear molecule is 
of course somewhat special, it is nevertheless interesting to study the dynamics of such a 
molecular system in the supercooled regime in order to check whether even such a simple 
system shows the phenomenon of a glass transition. Furthermore the continuous rotational 
symmetry around the long axis of the molecules simplifies also the theoretical description 
of this system, as compared to the molecules studied in Refs. p6 Hl8| , pO| -[2^ and thus allows 
to make a more stringent test of the theory. Although the MCT equations for a diatomic 
molecule in a simple hquid and for molecular liquids have recently been derived PB| , P7| , 
they have not been analyzed in great detail so far. Since the structure of these equations is 
not identical to that for simple liquids, it is not obvious whether the predictions derived by 
Gotze p8[ for MCT equations which apply to simple liquids are still valid. Thus one of the 
main motivation of the present work is to test whether these predictions hold also for the 
molecular system investigated here. 

The outline of our paper is as follows. In the next section we define the various correlators 
and diffusion constants that we will study. In addition we also summarize those results of 
MCT whose validity we will check for the present system. In section |T| the model as well 
as the details of the computer simulation are presented. The results are given in Sec. ^ 
and are discussed in the final section in which we also offer our main conclusions. 
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II. CORRELATION FUNCTIONS AND PREDICTIONS OF MCT 



In this section we define the correlation functions that are studied in this work and 
recapitulate some of the predictions of MCT. 

To study the static and dynamical properties of macroscopic systems with particles 
it is convenient to use correlation functions. For the translational degrees of freedom it is 
customary to characterize the dynamics with the help of the van Hove correlation function, or 
its space Fourier transform, the intermediate scattering function. For the rotational degrees 
of freedom convenient correlation functions are the functions Ci{t) and Cl'^\t) which are 
defined as follows: 

Cl{t) = ^T.(PliUt)-UniO))) 1>1 . (1) 



and its self part 



N , 

n,n' 



1 ^ 



n=l 



Here Pi{x) is the Legendre polynomial of order I and Un is the normalized orientational 
vector of the n-th molecule. These orientational correlation functions are straightforward 
generalizations of the intermediate scattering function to the case where also ODOF are 
present. More details on this point are given in Ref. |2^. The experimental relevance of the 



functions Ci{t) is given by the fact that for / = 1 and / = 2 they can be measured in dielectric 
and light scattering experiments, respectively. We also note that it is often assumed (cf., 
e.g., Ref. [^) that the cross terms in Ciit) can be neglected. In that case the experiments 
would also yield information on Ci^\ We will discuss this point below. 

Besides the intermediate scattering function or the correlators given in Eqs. (|l|) and @, 
the dynamics of the system can also be studied by means of the auto-correlation functions 
of the velocities Vnif) of the particles 

1 ^ 

m = J^T.^Vn{t) ■ (3) 

n=l 

or the corresponding one of the angular velocities oJn{t) 

1 ^ 

^Kt) = T7E(^Kcosa„(t))) />1 , (4) 

n=l 

where is the angle between uJnif) and ti}„(0). It was recently shown that ^2{t) is 

particularly useful to study the freezing of the ODOF |Q. Since the translational and 
rotational diffusion constants, D and Dr, respectively, can be obtained by means of a Green- 
Kubo relation from Eqs. (0) and (^), respectively, the significance of these time correlation 
functions is obvious. 

In the strongly supercooled regime and '^i{t) exhibit a negative, slowly decaying 
long time tail. This makes the numerical calculation of D and D^. from the Green-Kubo 
relation a difficult task, if the temperatures are low. Therefore it is customary to use the 
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corresponding Einstein relations, which are mathematically equivalent to the Green-Kubo 
identity. For the TDOF D is then determined from the mean squared displacement: 

1 ^ 

^= lim^E(l^n(t)-Xn(0)P) . (5) 

To obtain the analogon to this equation for the ODOF we replace Xnif) by the corresponding 
variable 4>n{t), which is defined as 

Mt) - 4(0) = A4(t) = /* dt'unit') . (6) 

Jo 

In analogy to Eq. (^) we thus obtain the following Einstein relation for the ODOF: 

1 ^ - 

Dr=]im-—J2{\Mt)-M0)\') ■ (7) 

Note that 0n(^) is not bounded, in contrast to Un{t), which is restricted to the surface 
of a unit sphere. This is the reason, why a replacement of Xn{t) by in Eq. (|^) would 

yield Dr = 0. 

Let us now recapitulate those predictions of the idealized version of MCT which are 
relevant for the present work. More detail can be found in Refs. The theory predicts 

that, in the vicinity of the critical temperature T^, all time correlation functions (j){t) which 
couple to the density correlation function should show a two-step relaxation behavior, i.e. 
exhibit a plateau-like region when plotted versus the logarithm of time. The time window in 
which (f) is in the vicinity of this plateau is called the /3-relaxation regime. The time window 
in which the correlator falls below this plateau is called the a-relaxation regime. 

MCT predicts that upon approaching Tc from above, the a-relaxation time r(T) diverges 
with a power-law, i.e. 

T{T)<x{T-ny^ , (8) 

with a critical exponent 7 > 1.5. Note that the values of T^. and 7 are predicted to be 
independent of the correlator. Furthermore the theory predicts that also the translational 
diffusion constant D shows in the vicinity of Tc a power-law behavior, i.e. 

D{T)(x{T-T,y , (9) 

with the same 7 as in Eq. (^. 

Finally the theory makes the prediction that the correlators should obey the so-called 
time-temperature superposition principle (TTSP). This means that if a correlator (j){t,T) is 
plotted versus the reduced time t/riT), the curves corresponding to the different tempera- 
tures fall, in the a-relaxation regime, on a master curve 0, i.e. 

0(t,T) = 0(t/r(T)) , (10) 

the shape of which is approximated well by a Kohlrausch- Williams- Watts function, i.e. 
0(t/r) ^ 74exp(— (t/r)^), where the amplitude A and the exponent (3 are not universal, i.e. 
will depend on (p. 
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III. MODEL AND DETAILS OF THE SIMULATION 



The model we investigate is a one-component system of rigid, diatomic molecules. Each 
molecule is composed of two different Lennard- Jones particles, in the following denoted by 
A and B, which are separated by a distance d and each of which has the same mass m. The 
interaction between two molecules is given by the sum of the interaction between the four 
particles which is given by the Lennard- Jones potential Vap^r) = ^eai3{{ao,p/ry^ — {cJap/rY} 
were G {A^B}. The Lennard- Jones parameters are given by aAA = (^ab = 1-0, ctbb = 
0.95, eAA = ^AB = 1-0 and esB = 0.8 and were chosen such that the system did not show any 
sign of crystallization even at the lowest temperatures investigated here. In the following we 
will use reduced units and use gaa as the unit of length, eaa as the unit of energy (setting 
kB = 1) and {a\xm/ ASeAAY^'^ the unit of time. If the atoms are argon-like this time unit 
corresponds to approximately 0.3ps. 

The choice of the intramolecular distance d between the A and B particles requires some 
consideration. On the one hand d has to be large enough to allow for a sufficiently strong 
coupling between the translational and rotational degrees of freedom. On the other hand it 
has to be so small that firstly the formation of liquid crystalline structures is unlikely ||31 



and that secondly the energy barrier involved in the intersection of two molecules is so large, 
that at the temperatures and on the time scale of the simulation such a crossing does not 
occur. We found that a value of 0.5 is a good compromise. 

In order to make the simulation more realistic we did it at constant external pressure 
Pe3;t=l-0. For this we equilibrated the system in the {N,p,T) ensemble, using the algorithm 
proposed by Andersen |^ and by setting the mass of the piston to 0.05. The length of 
these equilibration runs always exceeded the typical relaxation time of the system at the 
temperature considered, which allows us to conclude that in the subsequent production runs 
we were investigating the equilibrium dynamics of the system. After having determined 
from this equilibration run the appropriate density of the system for the temperature of 
interest, we fixed the total density to the so obtained density and started a production run 
in the microcanonical ensemble using the rattle algorithm |3^. Note that it is advisable 



to make the production run in the microcanonical ensemble if one wants to investigate 
the dynamics of the system, since the algorithms used for constant pressure simulations 
introduce an artificial dynamics which might lead to unphysical results. The step size we 
chose was 0.01 for the higher and 0.016 or 0.02 for the lower temperatures. For runs shorter 
or equal than 1.4-10^ time units these step sizes are sufficiently small to allow to neglect 
the drift in the total energy during the runs. This is not the case for the long runs at the 
two lowest temperatures which had a length of 3.0-10^ and 4.0 ■ 10^ time units, respectively, 
(= 1.5 ■ 10'' and 2.0 ■ 10^ time steps). During these runs the value of the total energy was 
reseted periodically (about 30 times during the whole run) to its value at the start of the 
run by rescaling the velocities of all the particles appropriately. Since this interference with 
the dynamics is only very weak it can be expected that the final result will essentially be 
independent of it. 

The temperatures we investigated are T = 5.0, 3.0, 2.0, 1.4, 1.1, 0.85, 0.70, 0.632, 0.588, 
0.549, 0.520, 0.500, 0.489, and 0.477. The total number of molecules was 500 and in order 
to improve the statistics of the results we averaged at each temperature over at least eight 
independent runs. 
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IV. RESULTS 



Before we start to present the results on how the ODOF freezes, it is useful to investigate 
the dynamics of the TDOF, since this allows us to estimate the temperature range in which 
the system is supercooled. Therefore we computed from the mean square displacement of 
the center of the molecules the translational diffusion constant D. MCT predicts that, in 
the vicinity of the critical temperature T^, the diffusion constant will show a power-law [see 
Eq. (^]. Thus we fitted our data for D with such a law, using Tc as a fit parameter. In 
Fig. |l| we show the diffusion constant versus T — Tc in a double logarithmic plot. Also 
included is the fit with the power-law of Eq. (^). We recognize that this fit is very good for 
a surprisingly large range in T and D. In particular this range is significantly larger than 
the one found for the atomic Lennard- Jones system []T1|,|T3|. Since no analogous analysis was 



done for the molecular systems studied in Refs. ||16 Hlq , p(]|] , and the range in D explored in 
Ref. |^T|] was significantly smaller than the one considered here, it is at the moment not clear 
whether the fact that the power-law can be observed over such an extended range in D is a 
peculiarity of the present system or whether it is a general feature of molecular systems. 

From the figure we also note that at the two lowest temperatures the diffusion constant 
is significantly higher than it would be estimated from the power-law which fits the data well 
at higher temperatures. The reason for this are likely the so-called "hopping processes" 



the contributions to the relaxation that are not accounted for in the idealized version of 
the MCT. Thus it can be expected that for such low temperatures the predictions of this 
idealized version of the theory are no longer valid. 

The value of the critical temperature Tc is 0.475 ±0.005 which allows us to conclude that 
the TDOF of the system become very slow in the vicinity of this temperature. This is also 
corroborated by our investigation of the intermediate scattering function (coherent as well 
as incoherent) for which we found that their a-relaxation time shows also a power-law with 



a critical temperature at 0.475 p3[. The critical exponent 7 of the power-law for D is 2.20, 
which is in the range of values found for 7 in simple supercooled systems. In passing we 
mention that for this molecular system the dynamics of the TDOF is qualitatively similar 
to the one of simple liquids and that therefore the molecular character of the particles 
does not seem to affect the dynamics of the TDOF significantly. 

We now focus our attention to the ODOF. The first quantity we investigate is the corre- 
lation function Ci\t) which was defined in Eq. (0). In Fig. |^ we show this time correlation 
function for all temperatures investigated. From this figure we recognize that for high tem- 
peratures Ci^\t) decays quickly to zero. At intermediate temperatures it starts to show a 
weak shoulder which, on lowering the temperature further, becomes more pronounced. The 
time range for which this shoulder is observed coincides with the one in which a plateau 
is observed in the intermediate scattering function |23| and thus can be identified with the 
/3-relaxation regime. 

From this figure we also recognize that for intermediate and low temperatures the shape 
of the curves in the a-relaxation regime seems to be almost independent of temperature, i.e. 
that the so-called time temperature superposition principle (TTSP) holds [see Eq. In 
order to investigate this point closer we plot in Fig. the same curves versus the rescaled 
time t/T[^\T), where t[ [T) is the a-relaxation time. We define t[^\t) to be the time it 
takes the correlation function to decay to of its initial value. From this figure we recognize 
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that the TTSP does not hold very well, in that the slope of the curves at t/Ti^\T) = 1.0 
changes significantly even at low temperatures. Thus we conclude that for this type of 
correlation function of the ODOF the TTSP does not hold very well, which is in contrast 
with the behavior of the TDOF of simple liquids [ll4| , p!5| and of the present system ||2^ as 
well as for C['^\t) for the OTP model studied by Wahnstom and Lewis |T^. We also note 
that defining the a-relaxation time Ti'^\t) in a different way, namely by the time it takes 
the correlation function to decay to 10% of its initial value, does not change this conclusion, 
since the TTSP does not hold with this second definition of Ti^^ either. 

Next we investigate the time and temperature dependence of C2\t), see Eq. (||), which 
is shown in Fig. |^. From the upper panel of the figure we see that (t) decays qualitatively 
similar to Ci'\t) (Fig. |^a). For the former, however, the height of the shoulder is lower than 
the one in Ci^\t), which is reasonable, since to a first approximation this height is given 
by the value of the second Legendre polynomial evaluated at the height of the shoulder in 
Cf^ [see Eq. (|)]. Note that this height is a measure for the corresponding nonergodicity 
parameter ff = limj^oo Ci(t) at Tc. 

Since the area under the a-peak and under the remaining part of the spectrum x'l (i-e. 
the critical decay and the microscopic peak) are related to /f and to 1 — fi, respectively, 
our result suggests that the minimum between the two peaks is less pronounced for x'l than 
for provided that the width of the microscopic peak is about the same for / = 1 and 
/ = 2. This could be the reason why the detection of this minimum is so difficult in dielectric 
measurements, i.e. / = 1, whereas it was readily found in light scattering experiments. 

From Fig. |^ we recognize that for this correlation function the TTSP holds well for 
times larger than T2^\t), the a-relaxation time for 6*2*^ (t), but that for shorter times quite 
significant discrepancies are observed, as it was the case for Ci'\ Thus we come to the 

(s) (s) 

conclusion that the relaxation behavior of CJ and C2 are qualitatively different. 

Since the time dependence of C[^^ seems to differ from the one of 6*2*^ we have also 

(s) 

studied the one of Ci for / = 3, . . . , 6, which, for the lowest temperature investigated, are 
show in Fig. ^d. From this figure we recognize that i) the height of the plateau decreases 
with increasing /, which can be rationalized by the same reasoning given above, ii) that with 
increasing / the correlators seem to show more and more a logarithmic time dependence in 
the /3-relaxation regime, and iii) that the correlators for odd values of / decay faster than 
the ones for even values of I. This effect can be understood by taking into account that at 
low temperatures the molecules make jump-like orientational flips of 180°, described in more 
detail below, which lead to a relaxation in Ci^^ if / is odd, but do not affect the correlators 
with even values of /. Furthermore we have found that with increasing value of / the TTSP 
holds better and better ||35| thus showing that from a qualitative point of view the correlators 
for the ODOF become more similar to the ones of the TDOF. 

In addition to the self parts, Ci^^ we have also determined the time dependence of Ci(t), 
one of the collective correlation functions of the ODOF [see Eq. (|I|)]. This correlation 
function is shown in Fig. |^ for all temperatures investigated. Although the noise in the data 
is signiflcantly larger than the one in C[^\ as it is often the case for collective quantities, 
we can clearly recognize that the time dependence of Ci is qualitatively similar to the one 
of C['^ and that the TTSP, cf. Fig. |, does not seem to hold. Furthermore we note that, 
e.g., at the lowest temperature and for t = 10^ time units, Ci'\t) is about 25% larger than 
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Ci(t), which demonstrates that the cross terms in Ci(t) should not be neglected, at least 
not in the strongly supercooled regime. 

From Figs. ^ and ^ we recognize that with decreasing temperature the relaxation 
of the ODOF slows down dramatically. Thus it is interesting to investigate the temperature 
dependence of the relaxation times. Since we have found that the diffusion constant (Fig. |I]) 
as well as the a-relaxation times of the TDOF show a power-law dependence on temper- 
ature, with the same critical temperature Tc, we checked whether also the a-relaxation times 
Ti'^\t), T2^\t) and ti(T) can be fitted with such a power-law (with the same Tc = 0.475). 
That this is indeed possible for about two orders of magnitude in r, is demonstrated in Fig. |^ 
where we show these quantities versus T — in a double logarithmic plot. As it was the 
case for the diffusion constant D, the values of r for the two lowest temperatures deviate 
from the power-laws, since also here the relaxation is too fast, which is likely to be related 
to hopping processes. Thus we conclude that these processes affect the ODOF also. 

From Fig. |^ we also note that the fitted power-laws do not extend to such high tempera- 
tures as they did in the case of the diffusion constant. This is evidence that, similar to simple 
liquids, the presence of a large temperature range for which such a power-law is observed 
is rather the exception than the rule. The critical exponents 7 of the three power-laws [see 
Eq. (H)], are 1.66, 2.42 and 1.52 for t['^\ 72'^^ and ri. Since the critical exponent for the 
diffusion constant is 2.20, we thus find that the four critical exponents are all different from 
each other, which disagrees with the prediction of MCT for simple liquids. However, if we 
determine the critical exponents for Cj^^^ for / = 3, . . . , 6 we find the values 2.25, 2.78, 2.55 
and 2.80. These values have to be compared with the critical exponent for the TDOF, which 
is around 2.6 |2^, and thus quite close to the one of C/*"* for the larger values of /. Thus 



this is more evidence that the latter correlators behave qualitatively similar to the ones for 
the TDOF. 

We also mention that a power-law fit to r/'*'' and ri, with the critical temperature 
as a free parameter, leads to a slightly improved fit. The optimal values of the critical 
temperature were to within 2% equal to 0.475, the value of Tc from the diffusion constant. 
Thus we find that the ODOF, measured by Cj;'^^ and Ci, would indeed freeze very close to 
Tc = 0.475, if the hopping processes were absent. 

Furthermore we note that a fit of r}*'* and T2^^ with the popular Vogel-Fulcher law, 
Aexp{B/(T — To)), also works remarkably well. In particular we find that this functional 
form is able to fit the data well at all lower temperatures, i.e. also the data points at the two 
lowest temperatures, for which the power-law failed to fit the data. Thus we conclude that if 
seen as a mere fitting function the Vogel-Fulcher law gives the better fit than the power-law. 
However, the Vogel-Fulcher temperature Tq was determined to be 0.328 and 0.386 for Ti^^ 
and T2^\ respectively. Thus we find that the two temperatures differ by about 20%, hence 
indicating that according to the Vogel-Fulcher fits there is no unique temperature at which 
the system ceases to relax. Therefore this sort of fit is, from a physical point of view, less 
appealing. 

The a-relaxation times Ti'^\t) and ti(T) are analogous quantities to the a- relaxation 
time t(T) of the intermediate scattering function. Since in supercooled liquids the temper- 
ature dependence of r and of the diffusion constant can be different (see, e.g., Refs. | ]T3| , pT| ) 



it is interesting to investigate also the rotational diffusion constant Dr, and compare it with 
the temperature dependence of the relaxation times Ti^^ and r2**'*. 
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As already mentioned in Sec. ||, the calculation of is numerically difficult when one 
uses a Green-Kubo relation. Instead it is much simpler to compute Dj, from the Einstein 
relation given by Eq. (J^). In Fig. |^ we show the time dependence of the mean square 
displacement of the angles i.e. ((A0(t))^) = (|(/>n(^) — 0n(O)P), where (pnif) — 0n(O) is 
defined in Eq. (H). 

From this figure we recognize that, analogous to the mean square displacement p3| , 
((A0(t))^) shows at short times a power-law with exponent 2.0, which corresponds to the 
free rotational motion of the molecules. For high temperature this type of motion crosses over 
directly into a diffusional behavior, i.e. ((A(/)(t))^) is given by a power-law with exponent 
1.0. This is not the case for the low temperatures, where the time regimes of the free 
rotation and the one of the diffusive behavior are separated by a time regime where the 
increase of ((A0(t))^) is slower than diffusive. The time at which this subdiffusive behavior 
starts is essentially the same as the one in which also the mean squared displacement of 
the particles (MSD) starts to show a subdiffusive behavior [^. In contrast to this the 



time where ((A0^(t))^) starts to show the diffusive behavior is, at the lowest temperatures, 
significantly less (by about 1-2 decades) than the times where the MSD starts to show 
the diffusive behavior. Thus, despite the qualitative similarity of the time dependence of 
((A0(t))^) and the MSD there are some distinct differences between the two quantities and 
thus we conclude that the plateau-like region in ((A0(t))^) should not be identified with the 
/3-relaxation regime. We will come back to this point later. 

From the time dependence of ((A0^(t))^) it is simple to compute Dr [see Eq. (^]. Note 
that, because ((A0^(t))^) reaches its diffusive limit at shorter times than the MSD does, 
the rotational diffusion constant can be calculated reliably from a relatively short run, an 
observation of which we will make use of below. 

In Fig. 1^ we show the temperature dependence of Dr in an Arrhenius plot. In order 
to facilitate the comparison between the rotational and translational diffusion constant, we 
have included the latter in the figure as well. (Note that we have multiplied D by 15 in order 
to make D and to coincide at high temperatures. Also it should be remembered that D 
and Dr have different units.) We see that for temperatures less than 2.0 the data (diamonds) 
can be fitted well with an Arrhenius law (solid straight line). Furthermore we recognize that 
the temperature dependence of Dr follows the one of D for high temperatures but that 
when the system starts to become supercooled, the curve for D drops significantly below 
the one for Dr- Thus we find that at high enough temperatures the ODOF and the TDOF 
couple strongly enough to show the same temperature dependence of Dr and D, which is in 
agreement with the well-known Stokes-Einstein- and Stokes-Einstein-Debye-relations. For 
lower temperatures D shows the power-law discussed in Fig. ^ the reason for which are 
likely the mode coupling effects. In contrast to this, Dr shows an Arrhenius law from which 
we can conclude that the rotational motion of the molecule is an activated process. We will 
study this process in more detail below. 

We have also checked whether at low temperatures, i.e. 2.0 >T> 0.477 the temperature 
dependence of Dr can be fitted well by a power-law and found that this is indeed possible with 
a critical temperature around 0.38 (dashed line in Fig. |]). This temperature is significantly 
lower than the critical temperature Tc we found for the diffusion constant, the intermediate 
scattering function and the relaxation times of Ci'^^ and Ci, which was 0.475. In order to 
discriminate between the two functional forms we made use of the observation described 
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above that the rotational diffusion constant can be determined from a relatively short run 
(compared to the a-relaxation times of the TDOF), see Fig. |^. Thus we used the temperature 
dependence of the density (obtained from our equilibrated runs at temperatures T > 0.477) 
to estimate the volume of the system at T = 0.41. We then set up the volume of the system 
such that its density corresponded to this extrapolated value and quenched the system to 
T = 0.41. After allowing the system to relax for 2.0 ■ 10^ time units we started to measure 
the time dependence of ((Ac/))^) for three independent runs. Note that this time is clearly 
not sufficient to equilibrate the system with respect to its TDOF, but it should at least 
allow the system to get reasonably close to its equilibrium state at this temperature. The 
so obtained ((A0(t))^) is included in Fig. |^ as well (bottom curve). We see that even at this 
low temperature the diffusive rotational behavior can be observed on the time scale of our 
simulation. Thus we could extract the corresponding value of and have included it in 
Fig. 1^ as well. This data point lies reasonably close to an extrapolation for the previously 
fitted Arrhenius law and is completely off the curve for the power-law. (The fact that this 
point lies slightly above this Arrhenius law can be understood by taking into account that 
at this temperature the TDOF are not quite relaxed. Hence it can be expected that all 
the relaxation times are smaller than they would be in an equilibrated sample and that 
therefore the measured diffusion constant is too large |5^.) Thus we conclude that the 
rotational diffusion constant, as defined in Eq. (|^, follows an Arrhenius law even at very 
low temperatures and that it is very unlikely that its temperature dependence is given by a 
power-law. 

In Fig. P we have seen that the relaxation times of the orientational correlation functions 
show a power-law dependence on temperature and that the critical temperature Tc is very 
close to the one of the translational diffusion constant or the one of the intermediate scat- 
tering function. From Fig. ^ it is recognized, however, that the rotational diffusion constant 
Dr does not show any exceptional temperature dependence in the vicinity of Tc. At first 
view these two facts seem to contradict each other, but as we will show now, this is not the 
case at all. In order for the time correlation functions C{ and C2 to decay to zero it is 
necessary that the orientation of the molecules changes by an angle on the order of 180° 
and 90° in the case of Ci^^ and \ respectively. In order to undergo such a large change 
of orientation, a molecule has to wait until the cage formed by the surrounding molecules 
breaks up. The time for this to happen is related to the relaxation time of the translational 

(s) (s) 

degrees of freedom and thus we find that the relaxation times of Ci and C2 become very 
large when the temperature approaches T^. 

For the rotational diffusion constant the situation is different. As can be seen from 
Eqs. (^ and (0), remains finite as long as there is a possibility that |A0(t)p increases 
(linearly) with time. At low temperatures a molecule will not be able to make large changes 
of its orientation but a small librational (tumbling) motion is still possible (for the TDOF 
this corresponds to the rattling of the particles in their cage), which was nicely demonstrated 
by Renner et al. for a system of infinitely thin rods on a cubic lattice |^ . It is not hard 
to see that this librational motion gives rise to a diffusive movement of the 2;-component of 
the vector and hence to a finite value of (here the 2;-axis is defined by the molecular 
axis of the molecule at time zero). 

In order to see this effect clearer we have investigated the orientational dynamics of the 
molecules at T = 0.41 in more detail. At this low temperature the orientation of most of the 
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molecules stays for a long time in the vicinity of the direction it was at time t = 0. Thus we 
determined the mean orientation of the 2;-axis of each molecule by averaging its direction 
over a period of 4 ■ 10^ time units. Note that such a mean direction makes only sense if the 
orientation of the molecule does not change significantly. Hence we will restrict the following 
analysis to only those molecules for which the maximum deviation from this mean axis was 
less than 45°. In the considered time window this is the case for 74% of the molecules. 
Having determined the mean z-axis, we chose a x and y axis perpendicular to the z axis and 
computed the time dependence of {{Acpait))^) = {\4>ait) — (/)a(0)p), with a G {x,y,z}. The 
time dependence of these three functions are shown in Fig. ^. We see that {{A(f)z{t)Y) is 
indeed significantly larger than ((A0a,(t))^) and ((A0y(t))^. This is in accordance with the 
picture put forward above that the orientational diffusion of the molecules in the z direction 
is much larger than the one in the x and y direction. We also recognize that the latter ones 
are not completely suppressed which is likely to be due to the fact that the cage in which 
the molecule sits is still slowly changing with time. It is important to notice that similar 
arguments do not hold for the TDOF. The rattling of the center of mass within a cage is 
isotropic an the average. It is the direction of the molecular axis that breaks this isotropy 
on a "mesoscopic" time scale. 

The fact that at low temperatures the molecules perform for a long time a librational 
motion can also be demonstrated nicely by considering the auto-correlation function \I'2(t) 
of the angular velocity uj{t), see Eq. (^). As was illustrated by Renner et al, '^2(t) is 
expected to show a plateau with a height equal or less than 0.25 if the motion of the 



molecule is of a librational type |^0[. In Fig. we show the time dependence of \l/2 for all 
temperatures investigated. From this figure we recognize that the short time relaxation time 
of \E'2 decreases with decreasing temperature and that at low temperatures the correlation 
function shows indeed a plateau. 

A different way to study the orientational motion of the molecules is to investigate the 
time dependence of the distribution function G{6, t), which is defined analogously to the self 
part of the van Hove correlation function, i.e. 

1 ^ 

Gi9,t) = j;^J2{6[9-cos-\m-um)]) , (11) 

i=l 

where Ui{t) is the unit vector parallel to the axis of molecule i at time t. In Fig. |TT]we show 
G{9,t) for different temperatures. Note that G{6,t) is defined such that for long times it 
approaches 1.0 for all values of 6. From Fig. |TT]a we recognize that for T = 2.0 this function 
decays monotonically in 6 for all times. This changes when the temperature is decreased to 
T = 0.63, Fig. |TTlb, since then, e.g., the curve for t = 77.7 (bold dashed curve) shows a small 
additional peak at 180° which is separated from the main peak at 0° by a shallow minimum 
around 120°. This additional peak stems from molecules which underwent a rotation of 180°. 
This feature becomes much more pronounced when the temperature is decreased further to 
T = 0.477 (Fig. pT|c). The mentioned minimum now exists for a large time range before 
it starts to disappear on the time scale of the a-relaxation. At even lower temperature, 
T = 0.41 (Fig. |lT|d), the minimum does not show any sign to fill up at all on the time scale 
of our simulation. However, we see that the peak at 180° is still observable, which means 
that a significant fraction of the molecules made a flip of 180°. 

To study these 180° jumps on a more microscopic level, we have also investigated the 
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time dependence of the angle 9 of individual molecules. From such studies we found that 
at low temperatures, i.e. T = 0.41, the 2;-axis of the molecules stays for a long time in the 
vicinity of its orientation at t = and then undergoes relatively quickly a flip of 180° (see 
Fig. |l^ for three representative trajectories). The typical time for this flip is around 50 time 
units, but also faster (see, e.g., Fig.|T2|b) as well as slower transitions can be observed. This 
transition time is relatively long compared to the time scale of a (translational) vibration 
of a molecule in its cage which is on the order of two time units. Therefore we conclude 
that such a 180° flip is not a fast process in which the molecule overcomes one barrier in a 
quasi-ballistic way, but rather the sum of a quick succession of smaller jumps. Finally we 
mention that the molecules do not show these sort of little jumps neither before nor after 
they undergo a 180° flip which shows that these little jumps are associated with the flips. 



V. CONCLUSIONS 

The main motivation of this paper has been to investigate the dynamics of the orienta- 
tional degrees of freedom in a supercooled molecular liquid. This was done by means of a 
MD-simulation for a very simple molecular system, a liquid of diatomic, rigid molecules. 

The first question we addressed was how the translational and the orientational degrees of 
freedom (TDOF and ODOF, respectively) slow down if the temperature is decreased and the 
system becomes strongly supercooled. Since the mode coupling theory (MCT) is presently 
the only microscopic theory which predicts, in its idealized version, a glass transition at a 
temperature T^, we have checked the consistency of our results with the predictions of the 
theory. 

In the a-relaxation regime MCT predicts in the vicinity of Tc a power-law behavior for the 
temperature dependence of the corresponding relaxation times and the diffusion constants 
[cf. Eqs. (^) and (^)]. We find that the translational diffusion constant D as well as the 
relaxation times r for the orientational correlators C\^'^ and Ci can indeed be fitted by a 
power-law and that the transition temperature is 0.475 ± 0.005. This fit describes the 
data for D and r very well for about four and two decades, respectively. The corresponding 
exponents are not universal, but vary between 1.52 and 2.45. 

The next interesting result is that the rotational diffusion constant D^, as defined by 
Eq. (0), shows a significantly different temperature dependence than the quantities just dis- 
cussed. For high temperatures and D are proportional to each other, in accordance with 
the hydrodynamic Stokes-Einstein and Stokes-Einstein-Debye relations. However, below a 
temperature T* ^ 1.4, which is far above Tc, D is described well by the mentioned power- 
law, whereas shows an Arrhenius dependence. We find that T* is the temperature at 
which the cage-effect starts to become important, i.e. the system begins to be supercooled. 
This can be inferred from the fact that at T* (Ar^) as well as ((A0)^) start to show anoma- 
lous diffusion behavior at intermediate times. Thus we conclude that the breakdown of 
this aspect of the hydrodynamic equations and the onset of the cage-effect in supercooled 
liquids occur at the same temperature. We stress, however, that the different temperature 
dependence of D and below T* is not related to the similar phenomenon observed in 
experiments In these experiments the relaxation of the orientational vector M„(t) of 



the n-th molecule is studied. The measured quantity corresponds to the correlator C\^\t) 
[see Eq. (pf)]. In most theoretical approaches the rotational dynamics is described by a 
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Smoluchowski equation in which it is assumed that the angular velocities can be eliminated 
adiabatically. This crucial assumption, and the subsequent linearization, leads to an expo- 
nential relaxation with a a-relaxation time r/^'' which is proportional to 1/Dj.l{l + 1) |p9|| . 
Using this relationship, can then be deduced. Our results at lower temperatures are not 
consistent with this theoretical result, because i) Dr{T) is not proportional to [r/'^''(T)]~'^ 
and ii) Ti'^\t)/t2^\t) ^ 3. The reason for this is the nonexponential relaxation we have 
found for which the proportionality of Dr and (t/*'')^^ and the relation T^^V'^'i'*'' = 3 is not 
granted. 

(s) 

In the vicinity of Tc the data for D and deviate from the power-law observed at 
intermediate temperatures. (The same deviations are found for the a-relaxation times of the 
coherent and incoherent intermediate scattering function of the center of the molecules [^). 
This is likely due to the hopping processes which restore ergodicity even at low temperatures. 
If hopping of the center of mass of the molecules becomes important at lower temperatures, 
this should be seen in the self part of the van Hove correlation function Gs{r,t) (see, e.g., 
Ref. []T3[)- Surprisingly, even at the lowest temperature, T = 0.477, Gs{r,t) does not show 
any sign of a second (smaller) peak at a distance r ^ 1, which corresponds to the mean 
distance between two neighboring molecules [^]. That D{T) nevertheless deviates from a 



power-law at low T may be explained by the jump-like reorientations of the ODOF which 
we have identified in the distribution function G{6,t) as 180°-flips. These jumps, which 
were also observed in earlier MD-simulations [|T^,0, may lead to a local "melting" of the 
neighborhood of a molecule which just jumped, whereby allowing a translational diffusion 
without the molecule having to jump over the walls of its cage. The rotational diffusion 
measured by Dr is an activated process, at least for lower temperatures. Its T-dependence 
can be described by an Arrhenius law even for a temperature T = 0.41 which is far below 
Tc. This Arrhenius dependence is somewhat reminiscent to the temperature dependence of 
the Johari-Goldstein /3-peak in the dynamic susceptibility ^0|. Since we have shown that 



the Arrhenius dependence of is related to the librational motion of the molecules, one 
thus might speculate whether the dynamics leading to the /?-peak is indeed related to such 
librations. However, in order to decide this one would have to investigate the equilibrium 
dynamics of the systems at significantly lower temperatures than it is presently possible. 

Regarding the time dependence of the various correlators at low temperatures we have 
found a two-step relaxation process, with a strongly nonexponential behavior, for all of them. 
This is in agreement with the results in Refs. []TB|,|TB|,^ and also with MCT. However, in 
contrast to the prediction of MCT, the time-temperature superposition principle does not 
seem to work very well for the orientational correlators investigated here, although it does 
so for Ci'^^ with / > 3 and for the correlators of the TDOF. This shows that, despite the fact 
that the temperature dependence of the relaxation times of the TDOF and the one of the 
ODOF are very similar, the relaxation dynamics of the two types of correlation functions 
is, if / = 1 or / = 2, qualitatively different. Finally we mention that at low temperatures 
C[^\t) and Ci{t) differ from each other in the a-relaxation regime by about 25%, which 
demonstrates that the cross-terms in Ci{t) should not be neglected. 
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FIGURES 



FIG. 1. Self-diffusion constant D versus T — Tc- The solid line is a fit with a power-law with 
the exponent 2.20. The dashed line is a guide to the eye. 

FIG. 2. Time dependence of C^'^ [see Eq. (2)] for all temperatures investigated, a) versus time t 
b) versus rescaled time t/T[^\T), where r}*^ is the a-relaxation time. 

FIG. 3. Time dependence of C^^ [see Eq. (2)] for all temperatures investigated, a) versus time 
t b) versus rescaled time t/Ti^\T), where r^^*^ is the a-relaxation time. 

(s) 

FIG. 4. Time dependence of for I = 1,...,6 for T = 0.477, the lowest temperature 
investigated. 

FIG. 5. Time dependence of Ci [see Eq. (1)] for all temperatures investigated, a) versus time t 
b) versus rescaled time t/T{T), where r is the a-relaxation time. 

FIG. 6. a-relaxation time t[^^ (squares), T2^^ (diamonds) and ri (circles) versus temperature. 
Solid lines: fits with power-law. The dashed lines are guides to the eye. 

FIG. 7. Time dependence of ((A^(i))^) [see Eq. (6)] for all temperatures investigated. 

FIG. 8. Temperature dependence of the rotational diffusion constant (diamonds) and of the 
translational diffusion constant D (circles). D is multiplied by 15 so that the two curves coincide 
at high temperatures. The straight solid line is an Arrhenius behavior and the dashed line is a 
power-law with a critical temperature 0.38. The dotted lines are guides to the eye. 

FIG. 9. Time dependence of ((A^^)^), {{A^yf) and {{Acf)^)^) and their sum ((Ac^)^) 
[see Eq. (6)] for T = 0.41. See text for details. 

FIG. 10. Time dependence of the auto-correlation function *2(i) of the angular velocities 
[see Eq. (4)] for all temperatures investigated. 

FIG. 11. The function G{6,t) for different times and temperatures. Consecutive curves are 
spaced by about a factor of two in time. The first curves (bold lines) correspond to a time of 
approximately 0.64. The inset show the same curves on an expanded scale, a) T = 2.0, b) 
T = 0.63, c) T = 0.477, d) T = 0.41. 

FIG. 12. Time dependence of the angle 9 of an individual molecule in the time rage where the 
molecule makes a 180° flip. Three representative curves are shown. T = 0.41. 
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